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We perform lattice simulations in pure-SU(2) Yang-Mills theory to investigate how the infrared 
behavior of electric and magnetic gluon propagators in Landau gauge is affected by temperature. 
We consider the largest lattices to date, in an attempt to keep systematic errors under control. Elec- 
tric and magnetic screening masses are calculated through an Ansatz from the zero-temperature 
case, based on complex-conjugate poles for the momentum-space propagators. As recently re- 
ported in [jl|], we find good fits to the proposed form at all temperatures considered, with different 
ratios of real to imaginary part of the pole masses for the longitudinal (electric) and transverse 
(magnetic) propagators. The behavior of the magnetic propagator Dj (p) is in agreement with the 
dimensional -reduction picture, showing infrared suppression (with a turnover in momentum) and 
violation of spectral positivity at all nonzero temperatures considered. The longitudinal propaga- 
tor Di(p) appears to reach a plateau at small momenta and is subject to severe finite-Af, effects 
around the critical temperature T c . As a consequence, only lattices with temporal extent N t > 8 
seem to be free from systematic errors. After these errors are removed, the infrared-plateau value 
is considerably reduced around the transition and the sharp peak observed previously for this 
quantity at T c is no longer present. The resulting infrared behavior for Di(p) at T c is essentially 
the same as for 0.5 T c . An investigation of the temperature range between 0.571 and T c reveals that 
a less pronounced (finite) peak may occur at smaller temperatures, e.g. T « 0.9T C . 
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1. Introduction 

At zero temperature, Landau-gauge gluon and ghost propagators are believed to be closely 
related to the confinement mechanism in Yang-Mills theories, following the so-called Gribov- 
Zwanziger scenario [§, |J]. A key feature of this scenario, formulated for momentum-space prop- 
agators, is that the gluon propagator D(p) should be suppressed in the infrared limit. Such a 
suppression is associated with violation of spectral positivity, which is commonly regarded as an 
indication of gluon confinement. Lattice studies (see [Q] for a review) have confirmed the suppres- 
sion of D(p) in the infrared limit and have also observed violation of reflection positivity for the 
real-space gluon propagator [Q]. The infrared data for D(p) are well fitted by a Gribov-Stingl form 
(see e.g. [||]), which generalizes the form originally proposed by Gribov, based on a propagator 
with a pair of complex-conjugate poles. These poles can be associated with complex values for 
dynamically generated masses, a behavior in agreement with the massive (or decoupling) solution 
of Schwinger-Dyson equations The same behavior is obtained in the refined Gribov-Zwanziger 
framework [§, g]. (See also [ JTot ] for a very recent proposal in maximally Abelian gauge.) 

At high temperatures, on the other hand, one expects to observe Debye screening of the color 
charge, signaled by screening masses/lengths that can in principle be obtained from the gluon prop- 
agator [ full . More specifically, chromoelectric (resp. chromomagnetic) screening will be related 
to the longitudinal (resp. transverse) gluon propagator computed at momenta with null temporal 
component, i.e. with po = (soft modes). In particular, we expect the real-space longitudinal prop- 
agator to fall off exponentially at long distances, defining a (real) electric screening mass, which 
can be calculated perturbatively to leading order. Also, according to the 3d adjoint-Higgs picture 
for dimensional reduction, we expect the transverse propagator to show a confining behavior at fi- 



nite temperature, in association with a nontrivial magnetic mass (see e.g. [12]). We note that these 
propagators are gauge-dependent quantities, and the (perturbative) prediction that the propagator 
poles should be gauge-independent must be checked, by considering different gauges. 

Although the nonzero-r behavior described above has been verified for various gauges and 



established at high temperatures down to around twice the critical temperature T c [12, 13], it is not 
clear how a screening mass would show up around T c . In the following, we try to use the knowledge 
gained in the study of the zero-temperature case to define temperature-dependent masses for the 
region around and below the critical temperature. We review briefly the existing lattice results 
for this temperature range, present our fitting form for the infrared region, show our preliminary 
results for gluon propagators on large lattices for several values of the temperature and draw our 
conclusions. A more detailed analysis and additional data will be presented shortly elsewhere [|4j]. 

2. Gluon propagators around and below T c 

The behavior of Landau-gauge gluon and ghost propagators around the critical temperature 
T c has been investigated in [|TJ]. That study showed a stronger infrared suppression for the trans- 
verse propagator Dr(p) than for the longitudinal one D^{p), confirming the dimensional-reduction 



picture also at smaller temperatures. [We note here that a recent study [16] discusses whether this 
suppression is consistent with Dj(0) = and investigates Gribov-copy effects for the propagators.] 
It was also found that the ghost propagator is insensitive to the temperature. For the longitudinal 
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gluon propagator, a very interesting behavior was seen: the data approach a plateau (as a function 
of the momentum) in the infrared region and, as a function of temperature, this plateau shows a 
sharp peak around the critical temperature. The exact behavior around T c (e.g. whether the peak 
turns into a divergence at infinite volume) could not be determined, since relatively small lattices 
were used. All studies mentioned so far are for SU(2) gauge theory. The momentum-space expres- 
sions for the transverse and longitudinal gluon propagators Dj{p) and Di(p) can be found e.g. in 
101. 



More recently, in [17], further simulations around T c confirmed the above results, and lattice 
data for the gluon propagator were used to construct an order parameter for the chiral/deconfinement 
transition. More precisely, the authors use a much finer resolution around T c and consider the 
SU(2) and SU(3) cases. A check of their calculation is done for the electric screening mass, taken 
as Dl{0)~ 1 ! 2 and extracted from the data, where only the p = raw data point is used. The con- 
sidered lattice sizes are still moderate. [We also mention a very recent study of the SU (3) case, 
presented in [jig).] 

Of course, even if an exponential fit of the (real-space) longitudinal gluon propagator works 
at high temperature, implying that Dl(Q)~ 1 / 2 is proportional to the electric screening mass in this 
limit, it is not obvious that this should hold at T > T c . One should therefore consider more general 
fits. At T = 0, the momentum-space propagator is well fitted by a Gribov-Stingl form (see e.g. 
[|]]), allowing for complex-conjugate poles 

This expression corresponds to two poles, at masses m 2 = a ± ib, where m = niR + imj. The 
mass m thus depends only on a, b and not on the normalization C. The parameter 7] should be 1 
if the fitting form also describes the large-momenta region (from our infrared data we get 77 7^ 1). 
For consistency with the usual definition of electric screening mass, we expect to observe mj — > 
(b — > 0) for the longitudinal gluon propagator at high temperature. Clearly, if the propagator has the 
above form, then the screening mass defined by D L (0)- l l 2 = y/(a 2 + b 2 )/C mixes the complex 
and imaginary masses iur and mj and depends on the (a priori arbitrary) normalization C. 



3. Results 



We have considered the pure SU(2) case, with a standard Wilson action. For our runs we 
employ a cold start, performing a projection on positive Polyakov loop configurations. Also, gauge 
fixing is done using stochastic overrelaxation and the gluon dressing functions are normalized to 1 
at 2 GeV. We take /3 values in the scaling region and lattice sizes ranging from N s = 48 to 192 and 
from N t = 2 to 16 lattice points, respectively along the spatial and along the temporal directions. 

We note that we have improved our procedure for determining the physical temperature T. 
Instead of using the value of the lattice spacing a in physical units to obtain T = 1 /N t a, we evaluate 
the ratio T /T c for a given pair (j3, A 7 ,) by expressing T in terms of the lattice string tension a. More 
precisely, we consider the ratio T /T c = y/G^/y/o, where a is the (lattice) string tension evaluated 
at (/3, N t ) and yfc~ c is evaluated at the critical coupling j8 c for the same A^. In this way, one avoids 
the inconsistency of obtaining different values for the physical critical temperature T c for different 
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p (GeV) p (GeV) 

Figure 1: Longitudinal gluon propagator at and around T c , for various lattice sizes and values of j3. Values 
for the temperature, x N t , /3, lattice spacing a and spatial lattice size L (both in fm, in parentheses) are 
given in the plot labels. Lattices with N t < 8 and with N t > 8 are shown respectively in the left and right 
panel. 



N t 's. We have taken the f5 c values for the various N t 's from [19]. The string tension has been 
evaluated using the fit given in [20]. This leads to slightly different values of the ratio T /T c as 
compared with [p], |5fl. 

All our data have been fitted to a Gribov-Stingl behavior, as described in the previous section 
(see Eq. 1A). These fits are shown here in all plots of Du{p), whereas a detailed discussion of 
the associated masses niR, mj will be presented elsewhere [14]. We generally find good fits to the 
Gribov-Stingl form (including the full range of momenta), with nonzero real and imaginary parts 
of the pole masses in all cases. For the transverse propagator Dj(p), the masses m« and mj are of 
comparable size (around 0.6 and 0.4 GeV respectively). The same holds for Di{p), but in this case 
the relative size of the imaginary mass seems to decrease with increasing temperature. 

Our runs were initially planned under the assumption that a temporal extent N t = 4 might be 
sufficient to observe the infrared behavior of the propagators. (Our goal was, then, to increase N s 
significantly, to check for finite-size effects.) For this value of N t , the chosen j3 values: 2.2615, 
2.2872, 2.299, 2.3045, 2.313, 2.333, 2.5058 yield temperatures respectively of 0.92, 0.96, 1.00, 
1.02, 1.05, 1.12, 1.98 times the critical temperature T c . (See comment above about slight differ- 
ences in values ofT/T c for a given /3 in our newest analysis when compared with [[[], |5|].) As seen 
in Fig. [I], the assumption that N t = 4 might be enough is not verified for the longitudinal propagator 
around the critical temperature, especially in the case of larger A^ . Indeed, as N s is doubled from 
48 to 96 and then to 192, we see in Fig. [l] (left) that the infrared value of Dl(p) changes drastically, 
resulting in a qualitatively different curve at N s = 192, apparently with a turnover in momentum. 
(Also, in this case the real-space longitudinal propagator manifestly violates reflection positivity.) 
We took this as an indication that our choice of N t = 4 was not valid and therefore considered 
larger values of N t . Note that, due to the improved method for introducing physical units, the data 
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T c 48° X 4 2.299 (0.17, 8.2) fm 
T c 96 3 X 4 2.299 (0.17, 16.3) fm 
T c 192 3 X 4 2.299 (0.17, 32.6) fm 
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Figure 2: Transverse gluon propagator at and around T c , for various lattice sizes and values of j3 . Values for 
the temperature, xN t , /3, lattice spacing a and spatial lattice size L (both in fm, in parentheses) are given 
in the plot labels. Lattices with A^ r < 8 and with N t > 8 are shown respectively in the left and right panel. 



for N t <% (left panel of Fig. |TJ) lie exactly at T c , while the data at N t > 8 (right panel) are at 0.98 T c 
and 1.01 T c . 

As can be seen in Fig. [j] (right), we obtain in this way a different picture for the critical 
behavior of D^{p). Once we use values of N t that are large enough — i.e. > 8 for T < T c and 
Nt > 8 for T > T c — the curve stabilizes within statistical errors for four different combinations of 
parameters. In particular, this includes the two curves at fixed physical volume (the blue and the 
magenta curves). It is interesting to note that the N s effects at T c are significant for N t = 6 (with 
opposite sign with respect to the N t = 4 case) and are still present for N t = 8 (and maybe also for 
N t = 16) slightly below T c , but not immediately above T c . Note also that the curves corresponding 
to the smallest physical spatial size (i.e. around 4 fm), may show mild finite-physical-size effects. 

In Fig. H we show our data for the transverse propagator Dj{p) at the critical temperature. In 
this case, the finite-physical-size effects are more pronounced and, in particular, the lattices with 
the smallest physical spatial size (the red and green curves on the right) show qualitatively different 
behavior when compared to the other curves. On the other hand, Dr(p) does not seem to suffer 
from the same small-A', effects as Dl(p). Also, we see clearly the strong infrared suppression of 
the propagator, with a turnover at around 400 MeV. 

In summary, the transverse propagator Dj(p) shows significant finite-physical-size effects 
at T c , while the longitudinal propagator Dl(p) is subject to two sources of systematic errors for 
small N t : "pure" small-A^, effects (associated with discretization errors) and strong dependence on 
the spatial lattice size N s at fixed N t , when this value of N t is smaller than 16. The latter effect 
was observed only at T < T c , whereas the former is present in a wider region around T c . For all 
investigated values of the temperature, D^{p) seems to reach a plateau at small p. 

The plateau value drops significantly for T >T C and then shows a steady decrease. The behav- 
ior of Dl(p) and Dj{p) for temperatures above T c is shown in Fig. [| We see, again, that Dj{p) 
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Figure 3: Longitudinal (left) and transverse (right) gluon propagator for T > T c , Values for the temperature, 
Nf xN t , j3, lattice spacing a and spatial lattice size L (both in fm, in parentheses) are given in the plot labels. 



shows finite-physical-size effects for the smaller lattices (especially for spatial sizes below 4 fm, 
but also for 8 fm when compared with 15 fm). In the Di{p) data, on the contrary, there are no 
visible systematic effects for the considered lattices at these values of the temperature. 

We also did runs at low temperatures, namely at T = 0, 0.25r c and .5T C . We show combined 
Dl{p) and Dr(p) data for these runs in Fig. |[ We see that Di(p) increases as the temperature is 
switched on, while Dj{p) decreases slightly, showing a clear turnover point at around 350 MeV. 
(Note that the runs at 0.5 T c on 48 3 x 8 lattices for /} = 2.299, 2.301 are equivalent.) As pointed 
out before in [jl]], the infrared behavior of Di{p) remains unchanged (within errors) from 0.5r c to 
T c . This can be seen on the bottom right plot of Fig. |], where we show (for clarity) only lattices 
with large enough N t and the largest physical size. More precisely, we take j8 = 2.299 on a 96 3 x 8 
lattice, corresponding to 0.5 T c , and j3 = 2.515 on a 192 3 x 8 lattice, corresponding to 1.01 T c . 
We see that, whereas the behavior of Dj{p) is consistent with a steady monotonic decrease with 
temperature, the fact that D^{p) stays invariant might suggest a flat curve for the infrared-plateau 
value of the longitudinal propagator as a function of temperature below T c . 

To investigate the issue, we have performed runs at other values ot T < T c . We have considered 
several values of T /T c , and studied the dependence of the infrared-plateau value with T /T c . In Fig. 
IH we show data for Di(0) for all our runs on the left-hand side, and for the region around T c on 
the right. We group together results from runs using the same value of N t , and indicate them by the 
label "DL0_A r ,". The data points indicated with "sym" correspond to symmetric lattices, i.e. to the 
zero-temperature case. Note that results for different A^'s at fixed N t may not fall on top of each 
other, which gives us an indication of the systematic errors discussed above. These are especially 
serious for N t =4 around T c (red points). We see that, surprisingly, the maximum value of Dl(0) 
is not attained for T = T c — as might have appeared to be the case from the N t = 4 lattices only 
— and it does not describe a flat curve from 0.5 T c to T c , as could be expected by looking at the 
bottom right plot in Fig. |j. Rather, it seems to lie at about 0.9 T c . Also, it clearly corresponds to a 
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Figure 4: Longitudinal and transverse gluon propagators at T = (top left), T = 0.25T C (top right) and 
T = 0.5T C (bottom left). Curves for T = 0.5T C and 1.01 T c are shown together for comparison on the bottom 
right. Values for xJVr, /3, lattice spacing a and spatial lattice size L (both in fm, in parentheses) are given 
in the plot labels, with the exception of the bottom right plot, which is described in the text. 



finite peak, which does not turn into a divergence as N s is increased at fixed N t . 

Finally, we also looked at the real-space propagators. We find clear violation of reflection 
positivity for the transverse propagator at all temperatures. For the longitudinal propagator, posi- 
tivity violation is observed unequivocally only at zero temperature and for a few cases around the 
critical region, in association with the severe systematic errors discussed above. For all other cases, 
there is no violation within errors. Also, we always observe an oscillatory behavior, indicative of 
a complex-mass pole. Typical curves for the longitudinal and transverse propagators in real space 
are shown (for T = 0.25T C ) in Fig. [6[ 
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Figure 5: Infrared-plateau value for the longitudinal gluon propagator [estimated by Di(0)] as a function of 
the temperature for the full range of T /T c values (left) and for the region around T c (right). Data points from 
runs at the same value of N t are grouped together and indicated by the label "DL0_A^", where "sym" is used 
to indicate symmetric lattices (i.e. T = 0). 




z (fm) z (fm) 

Figure 6: Longitudinal (left) and transverse (right) gluon propagator in real space for T — Q.25T C . Values 
for A^ s 3 xN t , j3, lattice spacing a and spatial lattice size L (both in fm, in parentheses) are given in the plot 
labels. Note that the solid lines are not fits. 
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4. Conclusions 

The transverse gluon propagator Dj(p) shows infrared suppression and a turnover in momen- 
tum (in agreement with the dimensional-reduction picture) at all nonzero temperatures considered. 
Also, it exhibits violation of reflection positivity as a function of real-space coordinates in all cases 
studied. The longitudinal propagator D^{p), on the contrary, appears to reach a plateau at small 
momenta, and does not in general show violation of reflection positivity. We have obtained good 
fits of our data to a Gribov-Stingl form, with comparable real and imaginary parts of the pole 
masses, also in the longitudinal-propagator case. This is in contrast with an electric screening mass 
defined by the expression D^O) -1 / 2 , which moreover may contain significant finite-size effects. 

The data for D^{p) are subject to sizeable discretization errors around the critical temperature. 
In particular, a severe dependence on the aspect ratio N t /N s is seen at T < T c for the smaller fixed 
values of N t . As a result, only lattices with N t > & seem to be free from systematic errors. After 
these errors are removed, we see an infrared value about 50% smaller than before. As noted in [jl|], 
this might suggest that there is no jump in the infrared value of D^{p) as T — > T c from below, since 
the resulting infrared behavior at T c is essentially the same as at 0.5T C (see Fig. ||). [We note that 
all previous studies of D^{p) around T c had employed N t < 4.] An investigation of the temperature 
range between 0.5T C and T c shows, nevertheless, a (finite) maximum of the plateau value Dl{0) at 
around 0.9 T c (see Fig. g). 
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